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Abstract. In this paper we prove convergence results for the 
boundary layer homogenization problem for solutions of partial 
differential system with rapidly oscillating Dirichlet data. Our 
method is based on analysis of oscillatory integrals. In the uni- 
formly convex regime and smooth boundaries we prove pointwise 
as well as If convergence results. Namely we prove \u £ {x)—Uq{x)\ < 
Q £ .(d-i)/2_i_^ Vi e D, V k > - 1, and for 1 < p < oo we obtain 
\\u E - u \\w < C K e K , Vk e (0, l/2p), where u e and «o are solutions 
of respectively oscillating and homogenized Dirichlet problems, 
and d(x) is the distance of x from the boundary of D. In particular 
for p = 2 our result relates to the recent result of D. Gerarad-Varet 
and N. Masmoudi IIGM1 . 



1. Introduction and main results 

Homogenization in partial differential equation is a well studied 
topic, and with a major impact in applications, in particular in ma- 
terial sciences, where impurities of material tend to spread all over 
and affect both qualitative and quantitative analysis of the materials 
properties. In situations as such, one tries to make an approximate 
averaging of impure quantities, to reduce the cost of any numeri- 
cal computations. The averaging and homogenization technique is 
now a well developed tool, with an abundance of literatures on the 
topic. Most of these literatures, so far, take on the analysis of the 
bulk equation, that relates to the heterogeneity of the material; e.g. 
highly oscillating flux, turbulences, singularly perturbed equations. 
We refer to the book of Bensoussan-Lions-Papanicolaou [ BLP[ for 
background and overview. 

The problem of homogenizing the boundary data for Dirichlet 
problem has long been a challenging task, and there are very few 
and incomplete results in this direction. We refer to a recent paper 
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of D. Gerard- Varet and N. Masmoudi [GMj for an introduction as 



well as some background in the topic; see also IlLSB . The problem 
of homogenization of boundary data has shown to be an extremely 
hard problem, with deep connections to ergodic and number theory, 
and diophantine approximation! 

To explain this to some extent, let us (for simplicity) consider the 
problem of finding a harmonic function in a (bounded) domain D 
in M. d (d > 2), with an oscillating boundary data g(x/e), where g 
is a continuous 1-periodic function in M d . Such a solution can be 
expressed by Poisson representation 

(1) u E (x)= I P(x,y)g(y/e)do(y), 

JdD 

where P(x, y) = n{y)A{y)S/G{x, y), with n(y) the outward unit normal 
to dD at y, G the Greens function and A the matrix representing 
the operator; see ©-©• For the Laplacian/harmonic case A is the 
identity matrix. 

If, instead of harmonic functions, we consider solutions to general 
operators of divergence type, then a similar representation is possible 
through the Poisson kernel of that operator. In particular our method 
is heavily reliant on such a representation. To our best knowledge, 
this approach was first introduced in fLSt 

The study of limit behavior of the solution now reduces to study 
of the integral above. In IlLSl , the authors use approximation, along 
with simple compactness and covering to show that when the sur- 
face dD does not have flat portions of positive area, with rational 
normals, then the limit exists and equals \dD\g. The identification of 
the limiting boundary value is an easy task, once one can use simple 
foliation geometry, and whether the so-called scaling surfaces cover 
the torus. This identification, nevertheless, becomes quite compli- 
cated if the operator is also oscillating, in other words if the matrix A 
in 0-® is periodic and oscillating A(x/e). 

It seems still an open question what the limiting boundary term 
should be. We refer to flGMl for such an analysis. 

Once the limit problem is identified, then the question of conver- 
gence rate arises. In HGML the authors present a convergence rate 



of the system (0, as e tends to zero. For an oscillating equation and 
boundary data, they prove that there is a limit solution Uq and a limit 
boundary data <p (with a complicated representation) such that 

\\u E - Uolb(D) ^ C a e K , V k < (d - l)/(3d + 5). 
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It should be noted that as the dimension tends to infinity, the order 
of convergence tends to 1 /3. Based on some analysis, we believe this 
speed should be 1/2, but we leave it as a speculation rather than a 
conjecture. 

It is our intention in this paper, to analyze the speed of convergence 
for a (non-oscillating) elliptic system, with an oscillating boundary 
data. The method presented here, does not seem to extend directly 
to the case of oscillating equations, and needs modification for such 
an analysis. However, for oscillating boundary data it seems that our 
technique is quite effective and quite likely to give the optimal rate. 
We could not achieve this here, due to extreme technicalities. Nev- 
ertheless, we are able to obtain reasonably good pointwise estimates 
which implies U estimates of order 1 /2p. 

It is noteworthy that our approach, which is based on harmonic 
analysis, and is of more general nature and it might even be of interest 
to Harmonic analysts. Readers may also consult several outstanding 
sources for the theory of homogenization and the references therein: 

BAAl, flJKOj, BAH, Em ETJ, fM- 



To state the problem at hand, let us start with fixing some notations, 
and definitions. 

Let D be a bounded convex domain in ]R d (d > 2), and T be its 
boundary. In this paper we study the asymptotic behavior of solu- 
tions to the following problem: 

(2) |div(A(x)VM £ (x)) = 0, inD, 

\u £ (x) = g(x/e), on I, 

where e > is a small parameter and A(y) = (A a ^(y)), 1 < a, $ < d, 

1 < i, j < N is a R N ' xrf2 -valued function defined on R rf . Using the 
summation convention for repeated indices the operator is defined 
as 



(3) L := -div [A(x)V] = - ' 



'J y 'dxP 



dx a 

We also consider the corresponding homogenized problem, namely 



(4) 



{ div(A(x)Vu (x)) = 0, in D, 
\u Q {x) = g, on T, 



where g = f g(x)dx. Due to results of |LS| solutions to (0) converges 



to HI) if, e.g., the domain is strictly convex. 
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1.1. Assumptions. We make the following assumptions: 

i (Periodicity) The boundary value g is 1 -periodic: 

g (x + h) = g{x), Vx e R d , Vh e X d . 

ii (Ellipticity) There exists a constant c > such that 

iii (Convexity) We assume that T is uniformly convex hypersur- 
face with all its principal curvatures bounded away from 0. 
The constants then depend on the smallest principal curva- 
ture. 

iv (Smoothness) For the boundary value g we assume C 1 smooth- 
ness. Also we suppose that all elements of A and the surface 
T are sufficiently smooth. Here we do not aim to obtain the 
optimal smoothness, but rather focus on the method itself. 

Our main results are the following theorems. 

Theorem 1.1. (Pointwise estimate) Let D c R d be a bounded, smooth, and 
uniformly convex domain. For each k > d - 1 there exists a constant C K 
such that 

(5) \u E (x) - u (x)\ < C x e^ l2 -j— t Mx e D 

where d(x) is the distance of x from the boundary of D. The constant 
C K depends on the smallest principal curvature, and other operator related 
ingredients. 

Using this theorem one can obtain If -estimates for 1 < p < oo, as 
in the following theorem. 

Theorem 1.2. (LP estimate) Let D c R rf be a bounded, smooth, and uni- 
formly convex domain. For each 1 < p < oo, and each k < l/2p there exists 
a constant C K such that 

\\U e - Uq\\u{D) < C K £ K . 

Ideas in the proof: The main idea in the proof is quite simple. Using 
the Poisson representation as in [TJ in conjunction with a Fourier 
expansion of the oscillating period function g(x/e) one realizes that all 
we need to do is to consider this as the Fourier transform of a surface 
carried measure, with singular weight. The problem thus reduces 
to studying the behavior of this Fourier transform. Such problems 
have been extensively studied in the past, in connection to Hilbert 
transforms, and singular integrals, se e.g. llPSB . It is nevertheless, 
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not so apparent how one can apply the already existing results and 
methods to the particular case of our integral. But to stay fair, and 
not take a complete credit for the methodology presented here, we 
have to stress that the novelty of our paper is the suggestion of the 
approach rather than the technique itself. These techniques are fairly 
standard in Harmonic Analysis, and are based on careful estimates, 
rather than neat use of geometry, that is more usual in PDE. 

We believe that our method can pave the way for further deepening 
into the subject, at least for the divergence type operators. It surely 
has the potential of being applied to other areas such as dynamical 
system and even number theory. But we leave this discussions out 
here, as the ideas are still very vague, and focus on the main problem. 

2. Preliminaries 

2.1. Notations. For e > set g £ (x) := g(x/e). For ieR i and r > by 
B(x,r) := {y e R rf : \x — y\ < r) we denote an open ball with center x 
and radius r. By P(x, y) we set the vector-valued Poisson kernel for 
©• 

Here we recall some standard multiindex notations. For a multiin- 
dex a = {a.\, a^) e Z d , and for a point x = (x 1 , Xd) s R rf we denote 
x a := x" 1 ■ ... ■ x" d d , for m = (m 1 , ...,m d ) e Z d we set \m\ := \mi\ + ... + \m^\. 
Also by ||ra|| we denote the Euclidian norm of m to avoid ambiguity 
with its modulus. 

In the sequel C, Q, C 2/ ... will denote absolute constants which can 
vary in different formulas. 

2.2. Preliminary tools. 

Lemma 2.1. For each a e there exists a constant C a such that 
(6) |D«P(x, y)\ < C g |x _^ +|ahl , x e D, y e I. 

Proof. Let G(x, y) be the Green's function of operator L := div(AV) in 
D, and let n(y) be the unit exterior normal to T at y. Then the Poisson 
kernel P(x, y) is defined as 

P{x, y) = n(y)A(y)V y G{x, y), xeD,ye T, 

hence it is enough to prove that 

|D-G(x,y)|<C a |y _^ +w _ 2 , xeD,yer, 

where C a depends on sup norm of G, and the operator L. 



6 HAYK ALEKSANYAN, HENRIK SHAHGHOLIAN, AND PER SJOLIN 

Fix x 6 D, y € T, and set r = \x — y\. Consider translated and 
scaled domain D := \{D - x), and set U = (B 2 (0) \ B 1/2 (0)) n D. Since 

G(z) := r d ~ 2 G(x,rz + x) will be a solution of the scaled problem, and 
since the domain U and all coefficients of the operator are smooth, 
then by standard elliptic estimates (see | ADN1 ] - 1 ADN1 1 ) we will 
have 

\D a z G{z)\ < C a , zeD, 

for some constant C a depending on a and independent of r. From 
the latter we obtain 

C a > |r d - 2 r |a| D';G(x,rz + x)\ = r d - 2+M \D^G(x,y)\ f 

hence the claim. □ 

Lemma 2.2. There exists a constant c > depending on space dimension 
only, such that 

\™ 2a \ > cIMI 2 *, 

ai+...+aj=k 

where k e N, m e Z rf , and a e Z+. 
Proof. Consider the function 

a^+...+a d =k 

where x ± 0. Clearly it is continuous and positive on the unit sphere, 

hence c := min f(x) > 0. Also observe that f(tx) = fix) for each t > 0. 

||*||=l 

Now the statement can be read off. □ 
Lemma 2.3. Iff 6 C k (T d ) and ft e R, then £ ^\c m (f)\ 

< oo, provided 

k + fi>\d. 



Proof. First observe that using integration by parts we get c m (f) 



(2nim) 

k. 



]^c m (D a (f)), for each m e Z rf \{0} andeacha e with«i+...+a d < 



Since D a {f) e L 2 (T d ) we have 

(7) Yu Z i^w)i 2 < «>• 



mtO ai+...+a tl =k 

From Lemma IZ2l and Holder inequality we obtain 

( \ 1/2 

E \Cm(f)\ <c y \ C m(f)\ ,, y 



m£0 mj=0 



|mi| 2ai • ... ■ \m d \ 2ad 



ai+...+a t i=k 
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m*0 



M/)l 



L 



\m.\ 



2«i 



2a d 



\m\ 



-(*+/}) 



< 



\a-[+...+a l j=k 



c 



xl/2 



m 



-2(fc+j3) 



£ £ |mf|c m (/)| 2 Yj 

^m£0 «!+... +a rf =fc / V m *0 

By virtue of © the first factor in the above product is finite. The 
second factor is also finite due to condition k + f} > d/2 of the current 
lemma. Hence the conclusion of the lemma follows. □ 

Now we formulate and prove our main lemma. Let£ = (£1, G 
]R d (rf > 2) be fixed, = (£1, and p be small positive 



numbers. Assume that for a function i/; : R rf 1 
C°(B(0,&)), tK°) = 0, |D^(z')| < C«, and 



R we have i/> G 



det 



dz'dz 

1 U 



>c>0, 



for each a G Zf" 1 and z' G B(0, fc). 

Assume further u : R d_1 — > ]R satisfies u G C°°(B(0,ao))/ supp(w) c 
B(0, |fl ), and |D fl w(z')| < C a for each a G Z^" 1 and z' G B(0,a ). 

Lemma 2.4. Under the assumptions and notations above for (d - 1)- 
dimensional integral 



u{z')dz' 



3 = 

\z'\<ao 

one has 

\d\ < cipm-v-w. 

Proof. Set A(t/i, t] 2 ) = p^) with A = (|£f + p 2 ^) 1 ' 2 and n \ + n \ = 1. 
Clearly A?]i = £'andAT? 2 = p^. Note that since A > (p 2 \£'\ 2 +p 2 £ 2 d ) 1/2 = 
p|£|, it is enough to prove that 

(8) |3| < CA" (rf - 1)/2 . 

Denote F(z') := r\ x ■ z' + rj 2 ^xp(pz'). Clearly VF(z') = r\ x + r/ 2 JV^(pz'). 
Then we have 



a= J" e^^'yzodz'. 



|z'|<fl 

Case 1. Vi/;(0) = and d = 2. 

For \r] 2 \ < C\, where C\ is a small constant, we have |F( Z ')I > c 2 > 0. 
Here we can invoke integration by parts in d, and using the fact that 
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the derivatives of u are bounded get an estimate \3\ < CA M , where 
M > is large, and hence also © for d = 2. 

If \t] 2 \ > ci we get \F"(z')\ > c 2 > since F'(z') = ^"(pz') and 
ip"(0) + 0. We therefore can apply van der Corput's Lemma (see flS]], 
p. 334) and obtain the estimate |#| < CA _1/2 . 

Case 2. Vi/>(0) = and d > 2. 

We first assume that l^l < C\, where C\ > is a small constant. As 
in the case d = 2 we integrate by parts and obtain \d\ < CA~ M where 
M > is large, and hence © follows. 

It remains to study the case when |7] 2 | > C\. We set Q p (z') = 
-iVi/'Ooz'Xandget 

VF(2 ' ) = " 2 (| + ^ (2 ' ) ) = " 2 (l" Qp<2 ' ) )' 

hence z' is a critical point of F if and only if Q p (z') = ^. Observe, that 
since the Hessian matrix of ip is non zero at z' = it follows that the 
mapping z' i-> y = Q p (z') is one-to-one close to the origin for every 
small p > 0. 
We have F(z') = r\if{z') where 

7 /2 P 

V/(z') = ^ + -V^(pz') = 5i - Q p (z'), 

and D 2 /(z') = D 2 i/'(pz / ), where D 2 denotes the Hessian matrix. We 
have 



3 = j e lniA ^ f{z,) u{z')dz' , 



\z'\<a Q 



where supp(u) c B(0, \a$). Since the determinant of the Hessian 
matrix of ijj is non zero at z' = there exists a ball B(0, R) and a 
positive constant C\ such that for any x' , z' £ B(0, R) one has 



1 

Ci' 



(9) - \x' - z'| < |Q p (x') - Q p (z')| < C a |x' - z'|. 



Clearly we may also assume that |D 2 /(x')| > c > for x' e B(0,R). It 
follows from ® that 

^|x' - z'| < |V/(x') - V/(z')| < d|x' - z'|, 
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and using the fact that Vi/^(0) = we obtain 

i |z'| < \Q p (z')\ < Qlz'l, 

for each x',z' G B(0,R). One can see that there exists a neighbor- 
hood M of such that if x' G M then there exists z' G B(0, |R) with 
Qp(z') = x'. Here the constants R, C\, and the neighborhood M are 
independent of r\ and p. Then choose a so that B(0,flo) c B(0, |R) 
andB(0,2Cifl ) C M. 

> 2da . We have 



First assume that 

|V/(z')| = 



^-Q P (z') 



> 



\Q P (z')\ > 



2Ciflo — Ci|z'| > 2Ciflo — = Ci^o/ |z'| < flo- 
Hence we can integrate by parts in 3 and obtain the inequality 

We then assume that 



< 2Ciflo- In this case — G M and there exists 
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z' Q G B(0, jR) such that Q p (z' ) = ^, that is z^ is a critical point for /. 
We have V/(Zg) = and therefore 

^ \z' - z' Q \ < |V/(z')| < C 2 |z' - z^| for z' G B(0,R). 

Now we can use Theorem 7.7.5 from 1 1 i] to obtain the estimate ©, 
and thus completing the proof when ip(0) = Vi/>(0) = 0. 
Case 3. Vi/;(0) + 0: 

mthiscaseweseti/>i(z') = i/'(z')-Vi/'(0)-z' / sothati/'i(0) = Vi/^i(0) = 
0. Further, for 

H(z') = ? -z' + ^-ipipz'). 

P 

we have 

H(z') = C-z' + £ d -0Mpz') + pV^(O) • z') = 
e ■ z' + ^V^(O) ■ z' + ^-i^i(pz') = + ^(0)) ■ z' + ^-^i(pz'). 



Next setting 



or 



V = 4' + W(0), 
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with c\E,\ < \v\ < C\B\, we shall obtain 

H(z') = v' -z' + v d -ip 1 (pz'). 

P 

Since Vi^i(0) = we arrive at 

\3\ < C(p\v\)-^ /2 < cipwrv-w, 

which completes the proof of the Lemma. □ 

3. Proofs of the main theorems 

Proof of Theorem 11.11 Take r > and consider a covering of 
T by the following family of balls 23 = {B(z, \r) : z e T}. From 
covering lemma of Vitali there exists a subfamily of disjoint balls 

M 

T> = \B{z, \r) : ] = 1,2, ...,M} c 2 so that T c JJ B{z jr r). Now fix a 

non negative (p £ C™(R d ), such that supp(<p) c B(0,2), and cp(x) = 1 
for \x\ < 1. For r > set <p r (x) = (p(x/r), so that supp((p r ) C B(0,2r), 
and for ; = 1,2, ...,M set <p r ,j{x) = <p r (x ~ Zj), so su PP(^r,/) c B(Zj,2r). 

(m y 1 

For ; = 1,2, ...,M denote B ; := B(Zj,2r) and <p ; - = I E 9r,nWj (pr,j(x). 
It is clear that each function (p; is defined on a neighborhood of T, 

M 

supp(<py) c B ; , and also £ <p ; (x) = 1, for xeT. Observe that from the 

definition of <py follows that for each aeZ 1 * there exists a constant C a 
depending on a and the function cp such that for all / = 1, 2, ...,M one 
has 

(10) \D a cpj(x)\ < Cr~ M , x e By. 

We also note that an easy volume argument shows that each point of 
T can be covered by at most 10 d different balls By. 

Now fix a small positive number a. Let x e D and consider a 
number < p < c rf(x), where c is a small constant, and d(x) = 
dist(x, T). We then take r = ap, so cpj has support in a ball By of radius 
lap. From ((10) we get 
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We have 

( M 



/IV 1 
p(x, y)[g E {y) -g\Yu WW 
U 1 



My) 



M r 

Yj I p foy)[&(y) -g\<Pi(y)My), 

7=1 n 



where T ; = T n B,. 

After a permutation of coordinates we may assume that there exists 
a constant b > and a smooth, real-valued function i/>(z') defined for 
\z' - Zq| < b, where z' = (zi, ...,Zd-\) such that 

r ; = {(z',^(z')): |z'-z(,|<10ap}nB ;> 

where B ; is a ball of radius lap. We may also assume that supp <pj C B ; > 
|D^(z')| < C fl , and \det{^)^{z')\ > c> for |z' - z^| < b. The last 
condition comes from the assumption of the surface being uniformly 
convex. 

Letting y^ e Tj we have 



I 



P(x, y)[ge(y) - g](pj(y)da(y) 



\x - y0)|« 



where 



r, = JV 



y<'»r i P(x,y)[ & (y)-J],>,(y)rf c r(y). 



It follows that 



I,- = J Ix-y^Y" 1 ^, (z', ^(z')))[g £ (z', ^(z'))(l+|V^(z')| 2 ) 

|z'-zj,|<10flp 

We now make a change of variable y' = z' — z' Q and obtain 

Ij = J \x - yUf^Pix, (z' + y', ip(z' + y')))[ge(z f + y' , + V')) ~ g\ 



\y'\<Wap 



<pfa + y'M4 + yO)(i + W^o + y'tf'^y' , 
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Wethenseti^i(y') = ip(z' + y')-So for \y'\ < b, where So = $(z' ). Hence 
i/>i(0) = and 

Ij = J \x- y^\ d - l P{x,{z' Q + y',s + Uv'))W + Y^o + W))-g\ 

\y'\<10ap 

<pfa + j/', so + iMy'))(i + NMy')\ 2 ) 1/2 dy' ■ 

We now set y' = pz' and obtain 

Ij = P d ~ l J \x- y^f-'Pix, (z' + pz', s + Mpz'Mgeiz'a + pz' , s + ^i(pz')) - 

|z'|<10fl 

cpj(z' + pz',s + i/>i(pz'))(l + W^ipz'ff^dz'. 

Set a = 10a. Since g is 1-periodic and differentiable we have g(y) = 
E c m e lnim -v. It follows that 

meZ d 

g £ (z' + pz', so + MP*')) = Yj c,^ 2 ™' [ ^- (z o + p z ' )+ ^ (So+ ^ ( p z ' ))] = 

meZ d 

y Cm g2n<[ x' z o+ ?s ]g2m[ § m'-z'+ ^ ( pz ')] 
meZ ri 

Hence 

ly = p'" 1 £ c,„e 2ra ^- 2 o + ? s ^ f |x- y^Y" 1 ?^ (zo + P z '' s o + <Mpz'))) 



n;#0 

|z'|<fl 



(py(z + pz',S + ^1(PZ'))(1 + IV^lCpz')! 2 ) 1 ^^"^'"'- 2 '^''' 1 ^'^^' = 

p^^c^ 2 ™^^ f <^ m '- 2 ' + ^<^»u(z')<fe', 

where 

«(zOHx-y Y~ 1 P(*,(z +P^So+^^^^ 

and m' = (m\, ...,m^\). It follows from the condition p < Cod(x) and 
Lemma 12.11 that for each a £ Z d + there exists a constant C a such that 
\D a u(z')\ <C a , for \z'\ <a . 
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Now by setting £ = jtn, = z m ' , and £d = ^ m d we obtain 



m*0 , ,r 



|z'|<«o 



where denotes the last integral in the sum above. From Lemma 
12 .41 we get 



m+Q mtO 

where the last sum converges by virtue of Lemma |Z3l We therefore 
get 



j P(x, y)[g E (y) - g](pj(y)do(y) 



<C 



1 



,(d-i)/2 



x - y0")| 



Hence 



M 



C£ ( d -l)/2 -0-1) y« Ejl < Ce (d-l)/2 -(d-l) f 1—doQ) < 

Ce (d-w -(d-i) r i i* - y\ 6 do{ ) < C£ (i-D/2 -ca-i) i 

r 

where 6 > is a small number. Here we used the fact, noted in the 
beginning of the proof, that for any y € T the number of / for which 
y e is bounded by 10 d . 
Now we take p = Cod(x) where Cq > is a small constant, and get 

1 



|tt e (x) - u (x)| < Ce (d ~ l)l1 

where 5 > is arbitrarily small. 
Theorem [Tj] is proved. 



d(x) d - 1+6 ' 



X G D, 



Proof of Theorem 11.21 First we consider the case p = 1. Using 
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Theorem [LI] we get 



e l/2 j 

<cjldt + cj e {d - l)ll t l - d - b dt< 



\\u E — WoIIl^d) 

£ l/2 

(11) Ce 1 ' 2 + Ce {d - 1)/2 e 1/2{2 - d - 5) = Ce V2 ~ 5/2 . 

Now assume 1 < p < oo. From ([TTjl and using the boundedness of u e 
and u we get 



\u e - WqIIlp(D) - 



^ \u e (x) - Uq(x)\\u e (x) - Uq(x)\ p x dx 



< 



c 



J \u e (x)- 



u (x)\dx 



<C(e 



,1/2-5/2 



flv 



C£ l/2p-6/2p_ 

Theorem II .21 is proved. □ 

Remark 3.1. By putting ip = 2 in Theorem 11.21 we obtain that for each 
k < 1/4 one has 

l|w £ - WoIIl 2 (d) < Ce K . 
Note, that the latter (for non-oscillating operators) strengthens the recent 
result of D. Gerarad-Varet and N. Masmoudi (see HGML Theorem 1) in 
dimensions up to 8, and in dimension 9 we get an equal convergence rate. 
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